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Abstract 

We review recent progresses on jet quenching measurements by PHENIX. With increased statistics, PHENIX has 
gone beyond the single hadron suppression Raa, and made measurements on multiple jet quenching observables, such 
as V2, /aa and Vj". We argue that, by combining these observables together, one can achieve a better understanding 
of the energy loss mechanism. We present new y-hadron correlation results with associated hadrons extended to 
low pt', an enhancement has been observed, suggesting a contribution of genuine medium response that is relatively 
unbiased by the initial geometry fluctuations. The status of full jet reconstruction and future perspective of PHENIX 
jet quenching program are discussed. 

Keywords: Quark gluon plasma. Jet quenching, single hadron suppression, di-hadron correlation, gamma-hadron 
correlation 



1. Introduction 

Jet quenching or suppression of high px jets and di-jets is one of the most intensely studied probes of the strongly- 
interaction quark gluon plasma (sQGP) created in Au+Au collisions at RHIC. Since the suppression level is directly 
caused by the interaction of the jets with the medium, jet quenching measurements allow us to infer the properties of 
the sQGP After nearly a decade long effort, jet quenching, as an experimental phenomena, has been firmly established 
at RHIC. However, the exact mechanism for jet quenching is still under intense debate. The often-used pertubative 
QCD (pQCD) framework, which works well in elementary p+p collisions, fail to describe simultaneously the light 
and heavy quark suppression. In contrast, non-perturbative approaches, for example those based on AdS/CFT gauge 
gravity duality |[T], seem to work well. Even within the pQCD framework, there are many models that are based 
on different and often uncontrolled approximations fT\, which, when tuned to the same data, predict very different 
medium properties. 

One way to improve the situation is to study multiple jet quenching observables at once, capitalizing on their 
different sensitivities to the energy loss mechanism. These observables include single hadron suppression R\\ and 
its azimuthal anisotropy relative to the reaction plane (RP) Raa(<P - ^rp) or V2, di-hadron suppression /aa and its 
azimuthal anisotropy relative to the RP /aa(0 - ^rp) or Vj", as well as modifications of electron-hadron correlation, 
y-hadron correlations and fully reconstructed jets. The idea is to fix the energy loss mechanism, while dramatically 
varying the path length that the probes transverse, for example by comparing Raa, V2, Iaa, and for light hadrons; or 
fix the path length and vary the interaction between jet and medium, for example by comparing di-hadron and electron- 
hadron correlations. The goal of this proceedings is to summarize recent PHENIX results on all these observables, and 
discuss the physics insights obtained by combining these results. Reader should refer to [3j for a theoretical review 
on this subject. 
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2. Leading hadron suppression: does Raa rise with 

Single hadron suppression Raa is the most studied observable for jet quenching. Due to the steeply falling spec- 
trum, this observable suffers from energy loss bias (observed hadron tends to have small energy loss), which limits its 
ability in distinguishing various model scenarios. Nevertheless, as the R^a measurements become more precise, one 
can regain some discriminating power. For example, recent precision measurements of ;7r" and i] mesons have started 
to pin down the shape of the /?aa at high p-p. This is important for understanding the energy loss mechanism. A slow 
rise of the Raa, for instance, is expected for radiative energy loss due to its logarithmic dependence on the initial jet 
energy. 

Results for the 0-5% centrality bin are shown in Fig. [l] where RUN4 jfi and RUN7 rj Raa El are fitted with a 
linear function Raa - b + m pr to extract the slope parameter m. The data do indicate a small positive slope, but the 
significance is less than Icr for RUN4 jfi and slightly above Icr for 77. The RUN7 tt" Raa seem to suggest somewhat 
more significant increasing trend with pr- The updated numbers should come out soon. 
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Figure 1: (Left) Raa for ff" 
and mesons from different 
runs in 0-5% centrality bin; 
(Right) 1 and 2 standard de- 
viation contours of tlie inter- 
cept vs slope for published 
RUN4 and RUN7 tj 
(The contour parameters for 
RUN7 Saa were con- 
verted to the function used 
for RUN4 via a simple lin- 
ear transformation). 



3. Azimuthal anisotropy of leading Iiadron suppression: wliat is tlie patli length dependence of energy loss? 

The overlap zone of the Au-hAu collision is not symmetric; consequentially, the emission rate of high pj- particle 
may vary with its angle relative to the reaction plane (RP) - *Prp. This azimuthal anisotropy can be characterized by 
either V2 parameter or RAAi<P ~ ^rv)^ whose magnitudes directly reflect the path length dependence of the energy loss. 

Figure|2]shows the pT dependence of the Raa separately for in-plane {R™p) and out-of -plane {R°aD- At low pj, the 
R'^ varies by almost a factor of 2 while that for out-of-plane direction is almost unchanged; this can be interpreted as 
a stronger radial flow influence in the in-plane direction. At high pj, the R°^j^ is more suppressed than R^aa^ reflecting 
greater path length for out-of-plane going jets. Interestingly, 7?°"^ also shows a stronger pj dependence than the in- 
plane direction, suggesting that the R°^ can better expose the true shape of the .Raa at high pT from jet quenching. 
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Figure 2: i?AA vs pj in 
the in-plane (left) and out-of- 
plane (right) directions, re- 
sults for three centrality bins 
are shown. 
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Figure shows the centraHty dependence of V2 at high pr \5\. The V2 values are large and increase toward 
peripheral collisions. They are compared to various pQCD model calculations \6]. These models are calculated with 
different geometry and different implementation of energy loss processes. For example, the HT and ASW models 
include only coherent radiative energy loss, while the AMY and WHDG models also include collisional energy loss. 
The calculations are tuned to the Raa value in the 0-5% centrality bin (right panel), thus the centrality dependence of 
V2 and Raa are predictions. All models describe the centrality dependence of Raa reasonably well, but significantly 
under-predict the V2 data. Furthermore, the calculated V2 values differ among themselves. Since Raa and V2 are anti- 
correlated, i.e. a small Raa implies a large V2 and vice versa, it is unlikely that one can describe both V2 and Raa by 
simply re-tuning the quenching parameters in these models. 




Figure 3: The centrality 
dependence of i'2 (left pan- 
els) and Raa (right panels) 
in 6-9 GeV/c range. They 
are compared with various 
pQCD models calculations 
(top panels), and schematic 
calculations that explore the 
roles of geometry and cu- 
bic path length dependence 
(right panels). 
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How to resolve this apparent discrepancies between data and theories? A recent estimation Q shows that the 
calculated V2 could be increased by 30-40% by considering modifications of collision geometry due to event-by-event 
fluctuations and CGC effects (dashed line in Fig. [3};). However, the calculation still falls below the data. Unless there 
are other modifications of collision geometry that we are not aware of, this large discrepancy implies that our current 
picture of pQCD-based energy loss is not complete. For example, the quadratic I dependence formula for radiative 
energy loss, AE oc dr^^pir, r + nr) oc drp^-^^iir + nr) ~ f (where ppart is participant density), may need to be 
changed to increase the anisotropy. There are three ideas along this line of reasoning. Liao and Shuryak [8] do this 
by requiring that most of the energy loss in sQGP be concentrated around Tc. One can increase the V2 by increasing 
the formation time tq to 1.5 - 2.0 fm/c |]9J. One can also increase the V2 with a cubic / dependence of energy loss: 

AE oc c/t ^2-^) p(t, r + nr) oc drrpp^yiir + nr) ~ P; such dependence is expected in certain non-perturbative 
energy loss calculations based on AdS/CFT gravity-gauge dual theory fTOl . Indeed, Fig.|3}; shows that our data agree 
well with calculations IfTTH Tl based on such I dependence. 

To further pinpoint the influence of the collision geometry and the path length dependence to high pr anisotropy, 
we studied the I scaling behavior of the Raa for various centrality and <p - *Prp angle bins. The idea is that if two 
selections have different A^pan and - 'I'rp, but similar average energy loss (AE), their suppression levels should be 
similar Thus if real energy loss scales as AE oc drr'"^^ pp^^ir + nr) = I„ (m — 1,2), then one expects Raa {hn)- 
The result of this exercise is shown in Fig.|4]for 6 centrality and 6 angles bins, that is 7?aa are plotted versus /„, at high 
Pt (7-8 GeV/c). I\{p^^^^ corresponds to quadratic I dependence in standard Glauber geometry, I\{p^^^^ corresponds to 
quadratic / dependence in CGC geometry with E-by-E fluctuations, and 7200^^^ corresponds to cubic / dependence 
in CGC geometry with E-by-E fluctuations. Comparing the three panels, we see that the Raa does not scale with 
hiptmt) ™^ ^i(PcGy' although the latter is slightly better; but it scales very well with hip^Q^^). 
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Figure 4: 

/,„ = I," dTT"'-'pp.Mir + nr) 
(m = 1,2) for 6 central- 
ity and 6 angle bins, cal- 
culated for standard Glauber 
geometry (left), CGC geom- 
etry with E-by-E fluctuation 
(middle and right). 
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Figure [5] shows the energy dependence of v'2 from ^/suN = 39 - 200 GeV. The V2 at ^/sf. 
gradual drop from 3 GeV/c to 7-10 GeV/c and remain positive at higher pr- The V2S at -^.snn = 39 and 62 GeV have 
limited pj reach, but their magnitudes are consistent with 200 GeV results, both at low pr <2 GeV/c where collective 
flow dominates and at pj- > 4 GeV/c where jet quenching should play a significant role. The former result suggests 
that, even at low energies, the medium already has to thermalize quickly and have small dissipation. The latter result 
is quite surprising, since we know that the Raa& at low energies are less suppressed than at -^snn = 200 GeV. It would 
be interesting to see how well the pQCD models, which describe the /?aa at low yfsMN lfT2l . can reproduce the V2- 
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Figure 5: The tt vi 
versus pj for yfl^N = 
200, 62 and 39 GeV. 



4. Di-hadron suppression I^a and its azimuthal anisotropy /aa(0 - ^rp)'- Further probing into the path length 
dependence of energy loss with away-side jet 

High px hard-scattered jets, at leading order, are produced in pairs that are back-to-back in the azimuthal direction. 
One observable of interest is the suppression pattern of the away-side jet opposite to the trigger hadron above a certain 
Pt threshold. This suppression is quantified by Iaa, the ratio of the per-trigger yield (away-side jet multiplicity nor- 
malized by number of triggers) in Au-i-Au coflisions to that in p-i-p colUsions. Pure geometrical consideration would 
suggest /aa < ^AA due to a longer path length traversed by the away-side jet. But a recent PHENIX measurement |fT3l 
shows that /aa is constant for associated hadron > 3 GeV/c, and this constant level is above the Raa for the trigger 
hadrons, i.e. /aa > ^aa (see Fig.|6]l. Furthermore, the constant level of /aa becomes even less suppressed for higher 
trigger pr- This result can be explained, at least partially, by the bias of the away-side jet energy by the trigger pr'. the 
initial away-side jet spectra are harder for higher trigger pr, consequentially, a larger fractional energy loss is required 
for the same /aa value. This result rules out the pure jet attenuation scenario where the jet survival probability depends 
only on the path length. 

The data in Figure |6] are compared with several pQCD model calculations. The quenching parameters of these 
calculations have been tuned to reproduce the Raa data, so the calculations can be regarded as predictions. The 
ACHNS model (based on ASW framework) tends to predict /aa ^ ^aa, thus disagreeing with the data, while ZOWW 
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Figure 6: Away-side per- 
trigger yield suppression /aa 
as function of associated 
liadron pr in central Au+Au 
collisions. 



model (based on HT framework) can describe the data rather well. This is another example showing that one can 
discriminate models by combing multiple experimental observables. 

PHENIX also measured the anisotropy of the away-side suppression, that is /aa as a function of angle relative to 
the RP (see Fig.|7]) 1 14 1. While the near-side /aa is essentially unmodified, the away-side /aa shows a strong variation 
with 05 = 0- ^Frp, sometimes by more than factor of two from in-plane to out-of-plane. This variation corresponds to 
an anisotropy parameter of Vj** = 0.29^q [j, and is much bigger than that for inclusive jt^ in the same trigger pr range 
v'2 = 0.13 + 0.01. This measurement is statistics limited, however if the result holds, it would have severe consequence 
for energy loss models since they usually predict much smaller anisotropy (for example, ASW calculation predicts 
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Figure 7: /aa as function 
of angle relative to reaction 
plane </>, = (()- in two as- 
sociated hadron pr bins (left: 
4-5 GeV/c, right: 5-7 GeV/c) 
at both the near- (circles) and 
the away-side (boxes). 



Before closing this section, we want to point out that the four observables, Raa, V2, /aa and Vj**, are intrinsically 
correlated. A smaller /?aa naturally implies a larger V2, and a smaller /aa implies a larger Vj**. The relation between 
Raa and /aa also depends on the away-side spectra shape controlled by trigger pj, but in general a smaller R^a 
implies a smaller /aa- Thus it is rather surprising to see that /aa is less suppressed than Raa (Iaa > Raa), yet has a 
larger anisotropy (Vj** > V2). 

5. Non-photonic electron-hadron correlation: probing tlie in-medium modifications of tieavy quark jets 

The observation of large suppression for non-photonic single electron (NPE), at a level similar to that for inclusive 
hadrons, remains a challenge for pQCD models ifTSl . Since most of these electrons come from semi-leptonic decay 
of charm and bottom mesons, the large suppression suggests that the charm and bottom quarks interact with the 
sQGP much more than expected from pQCD models. Non-photonic electron-hadron correlations provide valuable 
complementary information for heavy quark energy loss. This is because a large fraction of heavy quarks are produced 
in back-to-back pair at RHIC energy, and the fragmentation of the companion heavy quark is expected to contribute 
significantly to away-side hadrons associated with triggering electrons. 
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Results from a first measurement of the NPE-hadron correlation from PHENIX ||T6| (see Fig. [8]) suggests that 
the away-side hadrons are strongly modified relative to p+p. The away-side hadron yield, when integrated in a 
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Figure 8: Non-photonic 
electron-liadron corre- 
lation: (left) Per- trigger 
yield azimuthal angle 
distributions for p+p and 
Au+Au collisions, and 
(right) Iaa 



±30" window around A(p - n, indicates an enhancement below 1.5 GeV/c and a suppression above 1.5 GeV/c. 
Unfortunately, the large statistical uncertainties do not allow us to conclude whether the modification patterns are the 
same as that for di-hadron correlations fTTl . PHENIX has installed a silicon vertex detector (VTX) for the next run 
(RUNll). VTX has the ability to directly reconstruct and distinguish between charm and bottom mesons. This should 
allow us to directly correlate charm and bottom mesons with charged hadrons. We expect PHENIX to carry out the 
first measurement of charm and bottom separated heavy flavor suppression and correlation measurements in RUNll. 



6. Direct y-tiadron correlation: an unbiased probe for jet quenching and medium response 

Our discussions of the jet quenching so far have been focused on the single hadron and di-hadron correlation 
observables. These observables are subject to energy loss bias and geometrical bias, and they represent a complicated 
convolution of jets with different initial energy and different energy loss. In contrast, direct y and fully reconstructed 
jet are much less affected by these biases. For example, one can gauge the away-side jet energy with a direct y trigger, 
and one can systematically control the surface bias by studying the jet /?aa a function of jet cone size, y-hadron 
and jet-hadron correlations also give us direct access to jet modification and medium response. 

To leading order in pQCD, the energy of a direct y is a good approximation of the away-side jet energy. Thus one 
can measure the in-medium jet fragmentation via y-h correlation. Results for p-i-p and 0-20% AuH-Au 1 18| are shown 
in Fig. |9j plotted as a function of fragmentation variable ^ = - In {pjlp\- cos A^) x, - ln(z). Small ^ corresponds to 
high pj and vice versa. The fragmentation function from TASSO (mostly quark jets) and in-medium modified jet 
fragmentation function from a MLLA calculation, both for 7 GeV jets, are shown as lines to compare with the pH-p 
and Au+Au data, respectively. Note that the TASSO data and MLLA calculation are scaled down by factor of 10 to 
match the our data. This factor is needed since the PHENIX detector has limited rj acceptance, thus only catches a 
fraction of the fragments of the away-side jet. 

The /aa, or ratio of fragmentation functions of Au+Au to p+p, is shown as filled circles in the right panel of 
Fig.|9] The ^ range is cut off at 2 due to limited pj of the pH-p reference data. However, we can extend the ^ range 
by calculating /aa using the scaled TASSO data as reference instead (open circles). The two /aas are consistent with 
each other. The data clearly show a suppression at small ^ (large associated hadron pj) due to jet quenching, but an 
enhancement at large ^ (small associated hadron pj) possibly due to medium response to the quenched jet. 

The enhancement of the associated hadron yield at low pj was observed in di-hadron correlations [IT], as a dou- 
ble hump structure centered around one radian from n. However, the interpretation of this enhancement in terms of 
jet in-medium response, e.g. Mach cone, is complicated by possible contributions of E-by-E collective flow fluctu- 
ations [19|. It has been argued that such fluctuations lead to significant non-zero vi and V3 components which can 
mimic such double hump structure EOl . Since the direct y does not interact strongly with medium, it should have very 



small flow signal as indicated by PHENIX preliminary measurements on direct y v'2 (Left panel of Fig. 10 1. Hence the 
y-h correlation should be relatively unaffected by the v„ contribution, and the associated hadron at low pj is a robust 
measure for the medium response 11211 . Current y-h correlation indicates some broadening at the away-side (Right 
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panel of Fig. 10 1, however the statistical and systematic uncertainties are still too large for a definite conclusion. With 
the VTX installed in the next run, PHENIX expects to have a factor of 20 increase in the effective pair acceptance for 
associated hadrons around A<p = n/2,a region that is most crucial for medium response studies. 
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Figure 10: Left panel: the direct photon V2', 
right panel: direct y-hadron azimuthal dis- 
tribution compared with inclusive and decay 
y-hadron distributions. 



7. Modification of fully reconstructed jets 

PHENIX has carried out full jet reconstruction in p-t-p and CuH-Cu collisions using the Gaussian filter method Il22l . 
This method is infra-red and collinear safe, and is suitable for limited acceptance detector One of the primary 
challenges for full jet reconstruction in the heavy ion environment is how to handle the large underlying background 
fluctuations, especially at low pr- PHENIX employed a fake rejection method [22 J . in which a certain criterion is 
defined to suppress fake jets from background fluctuations that tend to have high multiplicity and large width. This 
method can directly suppress fake jets with high purity, so the systematic error due to underlying event is smaller than 
that for the direct background subtraction method. However, the jet sample passing the rejection criteria are subject 
to some efficiency loss and bias, which need to be evaluated carefully. 

Figure[TT|summarizes the current status of the jet reconstruction in PHENIX t22'|. The left panel shows the full jet 
spectrum in p+p collisions, unfolded to particle level, up to 60 GeV. The result is consistent with an NLO calculation 
and PYTHIA. The middle panel shows the jet Raa in central Cu-i-Cu collisions, unfolded to p-i-p jet energy scale, for 
two different jet cone sizes. Comparison between different jet cone sizes was argued to directly probe the jet shape 
modifications [23 J . We see that the R^a for larger cone size is less suppressed, but the uncertainty also is much larger 
presumably due to increased background fluctuation in a larger cone. The right panel shows the di-jet acoplanarity 
for several centrality classes. The widths extracted via Gaussian fits are consistent across all centrality bins, suggests 
a small kj broadening for surviving partons traversing the medium. 
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Figure 11: Left panel: jet spectra in p+p collisions. Middle panel: jet i?AA in central Cu+Cu collisions for two jet cones. Right panel: di-jet 
azimuthal correlation for several centrality selections in Cu+Cu. 



8. Future of jet quenching physics in PHENIX 

The primary goal for jet quenching physics is to obtain a coherent picture of the interaction of the jets with 
sQGP. The challenge for the field is that we are not yet able to simultaneously understand multiple jet quenching 
observables, such as ^aa^ Iaa, V2, Vj**, and heavy flavor suppression. In order to meet this challenge, we need 
not only more precise measurements for existing experimental observables, but also the capability to measure new 
observables that can provide much more detailed picture about jet medium interactions. Examples of latter category 
include reconstructed jets and di-jets in a broad acceptance and kinematic range, direct y-jet correlation, and heavy 
meson tagged jet, just to name a few. PHENIX has planned aggressive mid-term (20I0-20I5) and long-term (beyond 
2015) detectors upgrades to fulfill these requirements. 

In the mid-term, PHENIX will see the completion of VTX and FVTX detectors. These detectors should allow 
us to tag D and B meson directly, and provide extended acceptance for low pr charged hadrons for light/heavy 
hadron-hadron correlation measurements. The upgraded data acquisition (Super-DAQ) will take full advantage of the 
increased RHIC luminosity. In the long-term, PHENIX plans to replace the existing outer central detectors with a 
compact large acceptance EMCal and HCal, which together with VTX, FVTX, additional tracking layers and high 
DAQ rate, will allow us to measure jets, dijets, heavy flavor jets, and direct photon-jet correlations in a broad kinematic 
ranges. We refer more detailed discussion to li24J . 



9. Summary 

PHENIX has made several new measurements on single hadron and di-hadron correlation observables. By com- 
bining information from multiple observables, we are now able to better constrain jet quenching mechanisms and 
discriminate different models. We show that the 7?aa measurement, with increased precision, suggests a gradual in- 
crease at high pt- We find that the V2 at high pr exceeds the pQCD predictions, suggesting a non-trivial path length 
dependence of the energy loss. We also find that the /aa is less suppressed than Raa, Iaa > Raa, yet its anisotropy 
is larger than single hadrons, Vj** > V2- This result is rather non-trival given the anti-correlation between Raa and V2, 
and between /aa and Vj** . 

PHENIX also made good progresses on y-hadron correlations and full jet reconstruction. These measurements are 
challenging either due to their low rate (y-hadron) or large underlying-event background fluctuation (full jet recon- 
struction). By extending the y-hadron correlations to low associated pr, we observed strong evidence of enhancement 
due to energy dissipation of quenched jets. We have measured jet spectra in p-i-p and Cu-hCu collisions, and ex- 
plored the modification of jet shape and di-jet broadening. These measurements will benefit tremendously from future 
detector and luminosity upgrade of the PHENIX. 
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